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ABSTRACT 


Completion of the installation at the Naval Postgraduate 
Sencol of a cold-flow test facility for the turbine of the 
Space Shuttle Main Engine High Pressure Fuel Turbopump is 
፡ .ቲ ሰሪ] The article to pe tested is the first stage ot 
the "Alternate Turbopump Development' model designed and 
manufactured by Pratt & Whitney. The purpose of the 
facility is to enable the development of non-intrusive flow 
measurements and comparison of those measurements with 
numerical sinulations. Flow field characteristics of the 
Біле stator were predicted using a three-dimensional 
viscous flow code. A sensitivity study was conducted to 
determine the effect of inlet profile to flow field 
solution. Recommendations are made for future use of the 
test facility and validation of the numerical simulation 
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1. INTRODUCTION 


A. PURPOSE 

This thesis describes the modification, completion and 
initial test of a test facility designed to enable the 
development of non-intrusive measurement techniques within 
the rotating blade row of a transonic turbine. It also 
details the computer modeling of a turbine stator flow field 


Mara lidation of ап advanced viscous flow code. 


B. OVERVIEW 

The performance of a turbomachine depends heavily on the 
MiS fields generated within the machine during it's 
operation. Measurement, prediction and understanding of 
conditions in the flow field are essential to optimized 
design of turbomachines. More efficiently cooled blading, 
enabled by advanced computer modeling and design, can result 
in significant improvements in service life, reliability and 
lifecycle cost. 

Installation at the Naval Postgraduate School 
sKuwbobpropuls:on Laboratory (TPL) of the cold-flow test 
о (псу reported by Studevan [Ref. 1] was completed as part 
of the present work. The facility was designed to enable 
tests of a turbine intended to power the High-Pressure Fuel 


Turbopump (H”FTP) (Fig. 1) of the Space Shuttle Main Engine 
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Figure 1. 


(SSME). The turbine was designed and manufactured by Pratt 
& Whitney as the 'Alternate Turbopump Development' (ATD) 
model. 

The purpose of the new Turbine Test Rig (TTR) 
wenfiguration (Fig. 2) at TPL is to facilitate the 
advancement of non-intrusive Laser-Doppler Velocimetry (LDV) 
techniques required to validate viscous flow simulation 
codes. The long term goal of this research is to use LDV to 
measure the flow field immediately adjacent to the endwall 
of the turbine in the tip clearance region. An intermediate 
goal is to determine the performance of the single, and 
ultimately the two-stage turbine operating continuously in 
the TPL facility. Gaddis, Hudson and Johnson [Ref. 2] 
documented the exit conditions and overall performance of 
the ATD, while Hudson and others [Ref. 3] investigated the 
original turbine model. Both of these tests were conducted 
ша short daration blow-down facility at NASA Marshall 
Space Flight Center. 

1. Experiment 

Initial measurements of the ATD will determine the 
power output of the single stage turbine. The 
instrumentation documented herein will be primarily that 
concerned with the safe operation and control of the 
 ፡ ከ Су, Provisions have been made for initial flow field 


testing, which will be made with classical intrusive 
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Figure 


techniques. The initial tests were conducted to determine 
the integrity of the turbine and its subsystems. 

2. Numerical Simulation 

Numerical simulations of the flow field in the first 

stage stator of the SSME HPFTP turbine were accomplished 
using two codes developed by Roderick V. Chima called 
Turbomachinery C-Grid (TCGRID) (Ref. 4], and Rotor Viscous 
Ecce 3-D (RVO3D) [Ref. 5]. TCGRID was used for the 
generation of the 3D volume grid around a single blade 
element and RVC3D performed the numerical solution of the 


flowfield. 


II. FACILITIES 


A. EXISTING TEST CONFIGURATION 
The existing configuration of the TTR, as reported by 

Studevan [Reri. 1], eonsistedzer: 

* Inlet Ducting 

e ATD 

* Shaft and Bearing Housing 

+ Bearing Housing Suppext 

e Dynamometer 

* Bearing Lubrication System' 

* Controls and Instrumentation. 
All of the above listed components are contained within the 
TTR cell with the exception of the Controls and 
Instrumentation which are housed in a control room adjacent 
to the cell. The interior of the TTR cell may be seem from 
the control room through an explosionproof window. This 
allows the test rig to be operated and observed from a safe 
environment. The test arrangement is shown in Figure 2. 

During operation, high pressure air flows from the inlet 

ducting through the inlet strut housing and into the ATD. 
In the ATD, air passes through the annulus of the first 
Stage stator where it is accelerated and turned off axis. 


The air then impinges on the first stage rotor blades which 


move in reaction to the pressure, turning the rotor disk and 
converting the energy in the air to rotational energy. The 
air then exhausts into the TTR cell. As the rotor disk 
turns, it spins the main shaft in the bearing housing. The 
main shaft, via the quill shaft, is coupled to the 
dynamometer where the rotational energy is absorbed. 

Shaft speed of the turbine may be sensed and controlled 
by the dynamometer system. Configuration of the speed 
control system is unchanged from that reported by Kane [Ref. 
6]. Speed may be automatically or manually controlled by 
adjusting the water flow rate through the dynamometer. In 
the automatic mode, a magnetic speed-sensor signal is 
compared to a set point, and a controller drives the 
difference to zero. The dynamometer is instrumented to 
provide mechanical power output via torque measurements. 

The power oucput may be calculated in two other ways; 
namely, usinc the mass flow rate and temperature rise of the 
water through the dynamometer, or using the mass flow rate 
and stagnaticn temperature drop of the air through the 


то пе. 


B. CONFIGURATION CHANGES 
1. Bearings 
The original design of the shaft and bearing housing 
called for f»ur high-precision ball-bearing units. After 


coupling the shaft and bearing housing with the ATD, the 


assembly was to be dynamically balanced. Balancing proved 
to be impossible and the assembly was taken apart for 
inspection. Upon disassembly, one of the bearing units was 
found to be damaged. In order to simplify the system, the 
outer bearing spacer was redesigned so that only two ball 
bearing units were necessary (Fig. 3). Appendix A contains 
the engineering drawing. 

During reassembly of the shaft and bearing housing, 
the shaft was cooled in a refrigerator to 62'F and the 
bearings were slowly heated in an oil bath to 165 POEM 7 
noted that condensation formed on the shaft when it was 
removed from the refrigerator and it was wiped down prior to 
inserting the bearing. The shaft was placed on end on a 
workbench. When the bearing reached the proper temperature, 
it was removed from the oil bath with tongs and placed ona 
metal plate. From the plate, it was immediately picked up 
in gloved hands and carefully but quickly inserted over the 
shaft. When the bearing was inserted over the shaft in an 
attitude level with the workbench it did not bind and slid 
properly into position. The procedure was practiced several 
times with the damaged bearing before an attempt was made 
with the actual hardware. To insert the shaft and bearing 
assembly into the bearing sleeve, the sleeve was preheated 
on a hot plate to 100'F. The shaft and bearings were then 
pressed into the bearing sleeve with a milling machine used 


as a press. 
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Figure 3. Bearing and Shaft Showing Oil Mist Path 


2. Oil-Mist Cooling System 

The Portable Turbine Lubrication Unit described by 
Studevan [Ret. 1] for cooling the high-precision ball 
bearings was replaced by an oil-mist lubrication system. 
The system was powered by shop air which passed through a 5 
micron filter to remove particulate matter before entering 
the Norgren Micro-Fog Unit, where it forced a small stream 
of atomized oil through tubing to a nozzle (Ела. 4) ٦ 
nozzle was inserted into a hole atop the downstream end of 
the bearing sleeve (Fig. 3). The oil/air stream flowed 
through openings in the outer bearing spacer and impinged on 
the ball bearings. It then passed through the bearings and 
out two ports on the lower upstream and downstream ends of 
the bearing sleeve. This was an open system. The used oil 
vented to waste buckets and was not recirculated. 

3. Temperature Monitoring System 

Because inadequate cooling could lead to bearing 
degradation, temperature-sensing thermocouples were placed 
in holes through the bearing housing at each bearing outer 
race. These were wired to digital output meters installed 
in the control panel in the TTR control room and were used 
to constantl; monitor the bearing temperature and insure 
that the cooling system was operating effectively. The 
redesigned bearing housing assembly was rebalanced 


successfully at approximately 2500 rpm , with the oil-mist 
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Figure 4. 


1 


lubrication system operating. The temperature-monitoring 
system was installed, and the assembly was put into service. 
4. Vibration Monitoring System 

In order to ensure the safety of operating the TTR, 
a vibration-detection system was installed to detect 
possible out-of-balance conditions and bearing degradation. 
The vibration-sensing system consisted of four 
accelerometers and a signal analyzer. All four 
accelerometers were attached to the bearing-housing support. 
One was upstream at the top of the housing and the second 
was upstream but displaced 90° clockwise looking upstream. 
The remaining two accelerometers were installed on the two 
bearing-housing support struts. The first two were wired 
directly into a Scientific-Atlanta SD380 Signal Analy Zens 
the control coom. The two accelerometers installed on the 
Support strucs were wired into the control panel, and to an 
analog Ballantine Laboratories RMS indicator. From the 
analog indicator, they could be connected (if needed) to the 
Signal analyzer for digital processing and hard copy output. 

Prior to powering the turbine with air flow, it was 
analyzed to determine natural frequencies so that resonance 
at those frequencies could be avoided. The accelerometers 
were installed and wired to the signal analyzer before any 
of the inlet ducting was attached. The stationary turbine 


and bearing housing assembly was struck with a modally tuned 


L2 


hammer and the response of the accelerometers was recorded 
by the signal analyzer. An example from the top 
accelerometer may be seen in Figure 5. The top trace shows 
the real time response and the bottom shows the frequency 
spectrum with the first predominant mode at 67.5 Hz. The 
assembly was then turned by electric motor and a belt drive 
at varying speeds up to 960 rpm and the largest amplitude 
responses wete found and graphed, as shown in Figure 6. 
Trace A is from the top accelerometer and shows the 16 Hz 
speed signal along with its higher harmonics. Trace B is 
from the accslerometer at 90 , and shows similar 
characterist.cs. 
5. Control System 

The control system in the TTR control room included 
a Pusher-Porter electronic set-point controller for the 
outlet valve of the hydraulic dynamometer. This was used as 
a speed regulator with manual or automatic modes and was 
unchanged from the configuration described by Studevan [Ref. 
1] and Kane [Ref. 6]. In the automatic mode, the controller 
received a conditioned feedback signal from a magnetic speed 
sensor in the dynamometer. The controller compared the 
speed signal to a set-point and produced a corrective 
Signal. The corrective signal, received by the outlet 
valve, controlled the water level in the dynamometer. The 


speed varied in response to the resistance produced by the 
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Stationary Shaft Impact Response 


Figure 5. 
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Figure 6. 


water and the energy was dissipated into the water as heat. 
The manual inlet control valve was also located in the TTR 
control room. Manual operation of the speed control system 
was used during all the tests conducted in the present 
study. 

6. Probe Placement 

A probe access hole was cut through the outer casing 
of the turbine upstream of the first-stage stator. The hole 
location corresponded with the inlet of the grid designed to 
be used in the computer simulation of the flow field. A 
probe holder was designed and manufactured to hold the probe 
in place on the test rig. It was designed to be used in any 
position around the circumference of the turbine. However, 
a new access hole would need to be cut in order to use it in 
another position. An engineering drawing of the probe 
holder may be found in Appendix A. 

Combination probes (numbers 2 and 3) were placed in the 
inlet duct through holes cut in the piping. They were 
positioned at the same axial station to measure stagnation 
temperature and stagnation pressure into the turbine. 
Combination »robe number 4 was clamped to the table 
immediately downstream of the turbine to measure exit 


stagnation temperature and pressure. 
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III. EXPERIMENT 


A.  EXPERIMENTAL PROCEDURE 

Before any measurements were attempted in the TTR, the 
shroud inserts between the blade tips and the turbine outer 
casing were removed. With the shroud inserts removed, there 
was no danger of blade tips rubbing due to misalignment. 
Initial system validation runs documented herein were 
conducted with the shroud inserts removed. The TTR was 
operated twi-:e to gain experience, test safe operating 
methods and  shake-down' the facility. During the first 
shake-down run, speed control and vibration response were 
tested under various air and water flow settings. It was 
noted that if water flow was sufficiently reduced by closing 
the outlet valve to the dynamometer, the speed of the 
turbine climbed rapidly. On inspection of the TTR after the 
run, the dynamometer was found to be very hot. During 
inspection of the water-circulation system, water was found 
to have overflowed onto the floor of the pumping station due 
to surging d2veloped by changing demand to the dynamometer. 
It was decided that the water flow through the dynamometer 
should be set at the start of the run to keep the 
dynamometer cool and eliminate surging of the water 


Circulation system. 


TT 


For the second shake-down run, a specific procedure was 
adopted to set the water and air flow at the start of the 
run. In order to keep sufficient and consistent water flow 
through the dynamometer, the inlet valve to the dynamometer 
was set to 80% open and the outlet valve was set to 50% 
open. For mass flow to the TTR to be properly calculated, 
all the air rlowing through the metering orifice had to go 
to the turbine (Fig. 7). To ensure that all air flow went 
to the turbine, the number-two dump valve had to be 
completely ciosed. After the TTR inlet valve was set to 20% 
open, the nunber -two dump was closed, then the number -one 
dump valve was set to 40% open: Further manipulation of the 
mass flow to change the speed of the turbine involved 
opening the ITR inlet valve to increase speed or opening the 
number-one dump to decrease speed. It was necessary to have 
some air venting through dump valve number one in order to 
keep the total flow rate sufficiently high or the Allis- 
Chalmers compressor (Fig. 7) could surge due to excessive 
back pressure. 

During the experimental run, the procedure described 
above for setting water and air flow was followed. Pauses 
were made during the acceleration and system settings were 
noted. At each pause, the speed was allowed to stabilize 
and measurements were made of temperatures and pressures for 
the power oucput and mass flow calculations. Measurements 


relevant to power calculations are shown in Table I. 
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TABLE I. ABBREVIATED TTR DATA 


7 m5 ЭШЕ 





B. MASS FLOW AND POWER CALCULATIONS 

The turbine air mass-flow rate was metered by an orifice 
plate (Fig. 7) and was calculated using the method described 
by Vavra [Ref. 16] and summarized in Appendix B. The power 
output of the turbine was calculated using the stagnation 
temperature drop of the air across the turbine, and the mass 
flow rate. Table II shows calculated values of mass-flow 
rate, specific power out of the turbine, and efficience m 
Details of these calculations may be found in Appendix C and 


a complete table of measurements is contained in Appendix D. 
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TABLE II. TURBINE PERFORMANCE (SHROUDS REMOVED) 


Mass Flow Specific 


Rate | Power 
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IV. NUMERICAL SIMULATION 


A. GRID GENS RATION 

Numerical simulation of a flow field is highly dependent 
on the grid used to define the geometry, and its resolution. 
Many routines have been developed to generate grids within 
and about reulistic surfaces that bound fluid flows. The 
computer program used to develop the grids used in this 
investigation, Turbomachinery C-Grid (TCGRID), generates 
three-dimensional C- or H-type volume grids for 
turbomachinery. Grid files are compatible with RVC3D and 
PLOT3D [Ref. 7] which were used extensively in the flow 
56۱۲٢۴ 1٣6۲۰۰8 ы 

TCGRID generates a three-dimensional computational 
volume grid in the following manner: 


1. A course, equally-spaced meridional grid is generated 
between the supplied hub and tip. 


2. The blade geometry is interpolated onto the meridional 
SEDO. 


3. 2-D blade-to-blade grids are generated along the 
meridional grid lines in (m,rbar*theta) coordinates 
using an old version of the Sorenson GRAPE code 
[Бег а. 


4. The (m,rbar*theta) coordinates ате (стапвтоттеа не 
(2,1, theta сбоя Ее 


5. The 2-D grids are reclustered spanwise to make a full 
= G 
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— n | 


6. Finally, the (z,r,theta) coordinates are transformed 
Је cm О капут Ее ah. PLOTSD formar. 


Hub, tip and blade geometries must be in MERIDL [Ref. 9] 
۶07073٦ 

The hot-flow blade geometry data were received from NASA 
Lewis Research Center [Ref. 10] in (z,r,theta) coordinates 
appropriate tor use in TCGRID. Since the Naval Postgraduate 
School facility was cold-flow, a 99% thermal-shrink factor 
had to be applied to all linear dimensions. The FORTRAN 
code written to accomplish this may be seen in Appendix D. 

Grid generation and sizing is a subjective exercise in 
which there is no single 'right' answer. Often the success 
of a grid remains unknown until it is used with a flow 
solver. In this investigation, over one hundred flow 
solutions were attempted with various grid sizes before a 
successful combination was found. It was understood that 
the investigation required a fully three-dimensional viscous 
flow solution. Because of the resolution needed in the 
boundary layer, the grid had to be very fine. However , the 
successful grid contained more than twice the number of grid 
Peints in the spanwise direction than originally thought. 
The final grid size used for the flow solutions was 151 x 31 
x 45. Chima [Ref. 11] recommended 15 grid points in each 
endwall boundary layer for good resolution. An example of 


ehe TCGRID input file may be found in Appendix D. 
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Figures 8 through 11 show C-grids developed by TCGRID for 
the hub, midspan and tip. 

This geometry was difficult to model because of the 
endwall contour (Fig. 8). The hub and tip geometry was 
defined in tne input file with (2,г) coordinates. Then the 
stator contour was defined in five sets of (z,theta,r) 
coordinates. Each set of coordinates was at a constant 
radius. Since the endwall radius reduced at the stator tip, 
some of the coordinates were defined outside the endwall. 
TCGRID calculates the intersection of the stator coordinates 
with the endwall and generates the grid defined by the 
intersection contour. Figures 9 through 11 show grids 
developed for hub, mean and tip sections. Figure 12 shows a 
three-dimensional view of the final grid. The curved 
trailing edge shown in this figure also complicated the grid 


generation procedure. 


B. FLOW SOLUTION 
1. RUCIL 

Once a proposed grid has been developed, it can be 
used by a flow solver to simulate a flow.  RVC3D (Rotor 
Viscous Code 3-D) is a computer code developed at NASA Lewis 
Research Center for the analysis of fully three-dimensional 
viscous flows in turbomachinery. The code solves the thin- 
layer Navier-Stokes equations with an explicit finite- 


difference technique. It is applicable to annular blade 
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Figure 11. C-Grid at Stator Tip 
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Figure 12. Three-Dimensional Stator Grid 
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rows or linear cascades and includes two turbulence models 
along with a simple tip clearance model. Mathematical 
formulation of RVC3D is described in references 12 and 13. 

The code solves the Navier-Stokes equations in 
Cartesian coordinates with rotation about the x-axis. The 
equations are mapped to a general body-fitted coordinate 
system (Fig. 13). —Turbulence 15 modelled иза на аи та 
dimensional adaptation of the Baldwin-Lomax model (Ref. 14] 
or the Cebeci-Smith model [Ref. 15]. A sample RVC3D input 
file may be seen in Appendix D. 

2. Solution Method 
Simulations were run on two different computer 


systems; a Silicon Graphics Iris workstation and a Cray Y-MP 


EL-98 supercomputer. Both systems reside at the Naval 
Postgraduate School. Solutions were obtained for two inlet 
velocity profiles on each machine. The same version of 


RVC3D was loaded onto the two computers and identical input 
files were used. 

First, using the two computers for time efficiency, 
an inlet velocity profile to simulate a fully developed 
turbulent flow was found using the Iris while, concurrently, 
a profile simulating a slug flow inviscid profile was found 
using the Cray. The profiles were developed by modifying 
the RVC3D input files and running 100 iterations at a time 


to find the inlet velocity profile. After each run, the 
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Figure 13. General Body-Fitted Coordinate System 
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171) 


input file was modified and run again. The process was 
repeated until a satisfactory profile was obtained. It was 
assumed that the two computers would yield identical 
results. 

After several runs on each computer, it was noted 
that when similar values were used in the input files for 
the two machines, the velocity profiles that were calculated 
were vastly different. The turbulent velocity profile was 
run on the Iris until the solution converged and the same 
was done on the Cray with the slug profile. Then the input 
files were exchanged and the simulations were run again. 


All four solutions were plotted for comparison. 
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V. RESULTS AND DISCUSSION 


A. EXPERIMENT 

The data shown in Appendix C yield curves shown in 
Figures 14 and 15. The curves confirm general trends that 
were expected from the experiment. The data were taken 
sns S the Intention of producing a turbine map but to 
validate the data acquisition methodology and to check the 
integrity of the turbine and its subsystems. The data show 
that the syscem will work but that it is inconsistent. It 
was noted that the water manometers used to measure pressure 
drop across the metering orifice were contaminated. This 
explains the inconsistent values noted in Appendix C as P31, 
Е З апа P34. 

The difference in the stagnation temperature readings 
indicate that the flow may not have been circumferentially 
uniform. The stagnation pressure values PT2 are obviously 
erroneous and the data obtained with that combination probe 
were not used, 

The low efficiencies shown in Figure 15 were expected. 
Since the shroud inserts were removed for these runs, it was 
expected that much of the air flow would pass between the 


blade tips and the outer case. The air passing around the 
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rotor blades did not contribute to the power опери стии 
turbine and significantly reduced the efficiency. 

Figure 16 shows the speed sensor response of the TTR at 
800 rpm. The number in the lower left hand corner, measured 
in Hertz, comes from the rotational speed. But, since the 
magnetic speed sensor on the dynamometer detects the passage 
of teeth on a 30 tooth gear, the rpm of the system is twice 
the frequency shown. Figure 17 shows the speed sensor 


response at 7925 rpm. 


Б. COMPUTATIONAL FLUID DYNAMICS 

As noted previously, running RVC3D with the same input 
files on the two computers resulted in solutions with 
Clearly different inlet velocity profiles. Figure 18 shows 
the two turbulent inlet velocity profiles from one Mame 
file, and Figure 19 shows the slug flow profiles from the 
other. The differences in inlet profiles generated in the 
course of the solution, show that the inlet velocity profile 
is sensitive to either the architecture or the precision и 
the particular computer. 

Figures 20 and 21 show the differences in inlet veloción 
profiles using the two input files on the same computer. 
The similarity of the graphs shows that the output inlet 
velocity profile is not very sensitive to changes in input 
velocity profile. Comparison of Figures 18 and 19 with 


Figures 20 and 21 indicates that the output inlet velos 
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Figure 17. High Speed Vibration Response 
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Slug Inlet Velocity Profiles 


Figure 19. 
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Normalized Velocity 
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Pprorile 1S more Sensitive to computer than to input velocity 
profile. 

Flgures 22 and 23 show the midspan chordwise variation 
in normalized pressure distribution of solutions with the 
same input velocity profiles on different computers. 

Figures 24 and 25 show the same results with the different 
inlet profiles on the same machine. Comparison of these 
four graphs shows that all the variations gave identical 
pressure distributions. Therefore, the pressure 
distribution is sensitive to neither inlet velocity profile 
nor computer. 

Figures 26 through 29 are PLOT3D plots of midspan Mach 
number contours. The first two show the Cray solutions of 
the turbulent and slug flow inlet velocity profiles. The 
latter two show Iris workstation solutions for the same 
input files. Again, in spite of the differences in the 
input files and the velocity profiles, all four solutions 
are identical. The plots show the flow accelerating on both 
the suction and pressure surfaces of the stator with maximum 
Mach number in the area of the throat as expected. Boundary 
layer growth as well as wake shape and orientation are 
evident in the plots. 

It 1S important to note that the wake, in yellow, does 
Escorollow the contour of the grid but follows the exit 
angle of the stator. The cut that joins the pressure and 


Suction sides of the grid may be seen in Figure 10. The 
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Tact спас спе wake does, not follow this cut indicates that 
the solution is somewhat independent of grid geometry and 
ЕВЕ continu.ty is satisfied at each point in the flow. 

Figures 50 through 33 the total pressure contours across 
the exit plaie of the stator. These plots are also nearly 
identical bur. small variations do exist in the upper left 
corner of the exit plane. The turbulent solutions, Figures 
а 32, а-е different in this region from the slug flow 
сои опе, Figures 31 and 33. This would indicate that the 
solution is somewhat dependent on inlet velocity profile. 

In these figures, it is also possible to see the 
variation in wake location along the span. The yellow area 
at the edges of the plots shows that the wake bows toward 
the pressure side at midspan. This may be explained by the 
bowed shape of the trailing edge of the stator. As the flow 
turns around the stator, the extended arclength along the 
surface in the region of the midspan may turn the flow 
farther than the shorter arclength along the stator at the 
۲۶۰۰۰۰۰ 0.1 tip. 

In the boundary layer at the hub there is a small region 
of blue and yellow in the area of the grid cut. Although 
ПОШИ 1 scon:inuity, this locus coincident with the grid cut 
may indicate some solution dependence on grid geometry. 

An example of the residual histories of two solutions is 
shown in Figure 34. The plot shows that the residuals 


matched exactly for the first 3500 iterations. The Iris 
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„аъ 22723 == 


ТЕРЕНС | ۷ 0000 
HH BES ҮШІ 000/:0 
ИЛЛЕ НЕНИН ql JA 000CL 0 


رر گا لھگ ИТ И‏ تو سے وربور و ور 


WWW ON" m | | 009 


| MIA alió | 
| | ПАНА РАЗА 


| 
| 


| 
| 


| || 
A. 


07ج \ 


=== 


۱77113-11 


ҹә BN wet e па АМА 
— а 1 50411 


ا 04 90 3199 


m ud 

9ቱቱዐ[ቾሩ9'ፅ | Emp ውርር መርን i. eus | | 
я a. SS E a (ИГО 
OAC 000 — з )0()۱ 


2120 | ТСО Тее, 


Total Pressure = Cray - Slug Flow 


РИС Те) 
sourT 1nojuo?) o1nssoiq [£10], 50814 ]፤፤ ሻ ፤0)ዩ]5 IMAMI, dLAdH AWSS 


Figure ST 


ж — 








im 
Н 


== 
с 
mimi ut 


>» 


е ее е | am m Lm am um ае am uml a Й 


0 ЙИ | ши 
ШИН 0۳۳۵ 
ПИЛИ |۷۷ 
ا‎ ВИН НОМ 

| 00077۴۲ 


2 


= 


33338 
= oo Qs مر(‎ 
ете ላው 


Turbulent 


/ 00000 
| ШАР \ 
00۶۶۳ 


ч-ют 
= 


سے == 


== 

>) <> же) ج‎ 
<<< 
=== 


>> 
Г == 
: S 
ON e I stu 
CONONONONONONOO 


''ኸ 


|| M 
UAR C АА 
NAR UR M 


ris 


III 
ПАМ АН 


І 


ሙዴ --ኀ 
— 


* 


Со 
со 
O 
e 


\ 
WAY نت‎ ii ' 0 
SS И 


| AAN 
9 دن۹‎ 5 SAY 
== Thim 


= 
=" 
са 


ES 
cp X IE X IGT | u INN 
Є*+**01Х00`8 єс 
9»»01Х599 

Odd 00 0 x rH ኑትትትኑትሩ1 ее 


.+ 
9 ہر دیو دا ئ ج ۹ا 9 2 a‏ 


ZTZ0 EN 


Total Pressure 


| ۱7 آ لہ‎ 13117 11318107017, - vonnJoS Su] 
۲3017 10101770 ٥1015807۲ ]][07, 9001 117 1016018 3510717. 1ة‎ 7۶6 


Figure 32. 
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Figure 33. Total Pressure - Iris - Slug Flow 
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Figure 34. 


D 


Solution was continued beyond 3500 iterations in order to 
see the effect of varying the inlet velocity profile. Each 
of the variations of the input file is marked by a 
significant zise in residuals. However, after approximately 
500 iterations each solution converged to nearly the same 
residual value. Even after seven changes in the input file 
to vary the inlet velocity profile, the residuals converged 
to four orders of magnitude below their start. This shows 
that the solution with the input grid generated by TCGRID is 


very stable. 
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Ут. CONCLUSIONS AND RECOMMENDATIONS 


A. EXPERIMENT 

A working test model of the first stage of the SSME 
HPFTP ATD turbine has been put into operation in the Naval 
Postgraduate School's Turbine Test Rig (TTR). The TTR 
subsystems have been successfully exercised; however, all 
lines and connections should be inspected, cleaned and 
repaired as necessary before any measurements are taken for 
performance evaluation. Connections should be completed to 
the water temperature probes and torque meter on the 
dynamometer and data acquisition should be automated. Three 
methods of calculating the power out of the turbine should 
be automated and output on-line: namely, using the flow rate 
and stagnation temperature drop of the air through the 
turbine, using flow rate and temperature rise of the water 
in the dynamometer, and using the rotational speed and 
torque generated in the dynamometer. Once the system 
alignment is verified, the shroud inserts should be 
installed and the performance measured. 

The inle.- flow circumferential uniformity should be 
checked. If the flow is not uniform, a flow straightener 


should be designed and inserted into the inlet duct. It 1s 
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also recommended that the inlet and exit velocity profiles 
be determined experimentally 

Ultimately, the turbine outer casing will need to be 
modified to allow LDV measurements of the flow, particularly 


in the tip clearance region. 


B. NUMERICAL SIMULATION 

A working CFD model of the SSME HPFTP ATD turbine first 
stage stator has been generated. Four solutions were 
generated on two different computers. It is recommended 
that more work be done to investigate the inlet-velocity 
profile dependence on the computer system, as this is not 
acceptable. An investigation 72 needed to determine whether 
the anomaly at the grid cut on the hub is grid dependent, or 
is a real, physical phenomenon. 

In order to yield more realistic flow solutions, the 
inlet struts that hold the bullet nose upstream of the 
first-stage stator should be modeled so that their effect 
can be included in the flow solution.  Computations can be 
made with input of the actual inlet velocity profile, once 
this has bee: determined experimentally. For an analysis of 
the complete stage, the rotor geometry must yet be modeled. 
Such a modeling of the downstream portions of the turbine 
will allow comparison with experimental results when the 
turbine is tested with the rotor clearance at the design 


value. 
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APPENDIX B. CALCULATIONS 


A. MASS FLOW RATE 


The following assumes the 9.0 in. orifice is installed 


and is derived from Vavra [Ref. 16]. 


E 
نج‎ 265 
Ti 
Р 
ሶፎ, m 
Ti 


‚_ / 
ВЕ ОВ ДЕ Е 
Z 


Perse solve: 





with (=1; 





NOW папа: 


where, 





m=C th’ in lbm/s. 


t,= temperature ahead of orifice (°F) 
h,= pressure drop across orifice (in. H,O) 


P,= absolute pressure ahead of orifice (in. Hg) 
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1. Flange Taps 


177 


T, 





m'-ca Y, 


Where 7. 0.4002 


Therefore, 


t.-60 h Р 
ih/-6.6445 [140.0015 2289) т О, E 
100 Р, N Е, +460 


Е 
= 1. ስደ ሰ 5 ንሳ) 
100 





‚ / 
х=0.2525 1 
Z 


(214 0.0047 
X 
m=Çç m, 


2. Vena Contracta Taps 


Е,-60 h DES 
ео 0291 — m 
100 P, ር, +460 


С 
100 





ОНО 2 21)‏ و 


‚ / 
=e oe Ge 
Z 


‚ 0.0029 


1 
: x 


m=C rh’. 
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B. TURBINE POWER 


Power=mC_AT 


8 9 4 9٤9 Еси Ре: ur. 


AT=TT35= Tra 
TT4 = Stagnation Temperature of Air at Turbine Exit (`F) 
TT3 = Stagnation Temperature of Air at Turbine Entrance ('F) 


Since 1 Hpos = 550 ft.lbf and 778 ft.lbf = 1 BTU, 


Power=m" 0.24 (TT3-TT4) 222. 
550 
C TURBINE EFFICIENCY 
peu 
F joo) y O 
Y Pp =p 
where TT,, = Isentropic stagnation temperature is given by 
TISTI (yo 
4i 3 РТ 
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COMPLETE TABLE OF MEASUREMENTS 


APPENDIX C. 
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APPENDIX D.  FORTRAN PROGRAM AND INPUT FILES 


RVC3D Input rile 


"SSME HPFTP -- Indigo -- Real Slug Flow' 

6111 1п5151 jm=31 km=45 itl=20 iil=-70 &end 

&nl2 cfl-4.0 avisci-0.0 avisc2-0.5 avisc4-1.5 ivdt-1 nstg-4 
itmax-985 

irs-1 epi-0.50 epj-20.60 epk-0.60 &end 

&nl3 ibcin-3 ibcex-3 isymt-0 ires-1 icrnt-50 

iresti=1 iresto=1 ibcpw=0 1а1п=0 «епа 

&nl4 emxx=0.212 emty=0.0 emrz=0.0 expt=0.0 prat=0.6500 
q А 

om=0.000000 igeom=1 alex=-60.0 &end 

&nl5 ilt-2 tw-1.00 renr-6.651e6 prnr-.7 prtr-.9 
vispwr=.666666 

srtip-0.0 cmutm-14. jedge-15 kedge-11 iltin-0 dblh-0.00 
ab1t=0.00 £end 

&nl6 101=1 102=151 oar=0. ixjb=0 njo=1 nko=3 

jo=1 ko=13 17 21 «епа 


> = - = - - т - - - - — - - - - - - - - = = - - - = መ መ መ መ = መ መ መ  ፎ<- መ መ መጨ - ~ - መኗ መ <= መ መ ~ 一 መ መ = = መ መ = - - - -- 


6:7 


Inlet.f - Reads RVC3D output and builds files for plot Tukin 
Inlet Velocity Profiles. 


QU kk ck kk k  k  k X X k kk kx kk kk C kk kk kk پر پا پا پر پر پر پر پر پر پر پا پر پر پر پر پر پر پا‎ A پر‎ A پر پر پر پر پر بل‎ ЭХ 


= inlet.f reads rvc3d files & writes ascii files for 
рони 

се unit 1 = input xyz file Eo ۶۹۹ и 

С unit 3 = Input су ше (TOTE: 3) 

E unit 17 = output inlet velocity profiles 

E unit 4 = output residual history 


СХХХХАХАХЗХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХЭ 


parameter (ni=-152,njJ=32,nk=46) 
real x(ni,nj,nk) y (mi nj mk pc maw mh 
real ےو.‎ )5 , 37,673,۰۴۴۶ 6 ٦ 
real vel(5) z2 (5) 
dimension ii(5) 
E i-values are hard-wired below 
data 11/60, 70. 725. 69909. 


СХХАХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХХ a XK 


С read g%1d coordinates 
C Kk Kk k k k k Kk k k k k > k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k 


read (1,*)im,]m,km 


теаа(1,*) СООО Еа 

1 «(су (1,3.,К1, 261 , 10202, Та ٥٠ 

2 СС ЈК Е т ٣٣ى‎ 
QU CX C Ck Ck kk Ck kk Ck kk Ck Ck C CC CCCII CIC CIC C000 CIC Ck 0C Ck kk Ck CX Ck kk A ko ko X koX ko ok XK 
ር: read restart file 


KKEKKKKKKKKKKKEKKKKKKKKEKKKKKKKKKKKKKKEKKKKKKKKKK KK KKK KKK A I I И. 


read (3,*) imax, Jmax, kmax 
read(3,*)fsmach,alpha, re, time 


e 
icheck-iabs(im-imax)-*iabs(jm-jmax)-*iabs (km-kmax) 
if (icheck.ne.0)then 
write(6,610)im,jm,km,imax,Jjmax,kmax 
stop 
endif 
© 
read(3,') ((((qq(1,1,3,k),21=1,1m) ,j=1,3m) , k=1, kn) ae 
С 
С additional residual data 


read(3,*)itl,ill,phdeg,ga,om,nres,dum,dum,dum,dum 
read(3,*) ((resd(nr,1),nr=1,nres),1=1,5) 
C kk kk kk چاو پل‎ እና ደ እ ቷ KKKKKKKK KK AAA AA AAA AAA AAA AAA XK 


e vel to unit 17 

QC XC C C C C C Ck Ck Ck CC چاو چاو چاو‎ Ck Ck kk Ck kk kk kk Ck kk kk ko ko ko ck ko ko ko ko ko Ck ko k ko ko ko k k k X 
itr- iml itl 
ggm-ga*í(ga-1.) 
0 707 

С normalize x by chord 
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соо ив 
k=kk (1) 
зле жк ла рок 
Шах Е.К) 
do 5 а ЕЪЕЁ ТЕТ 
ев (ем, 7 ድሮ) 
Бе “шах та нах, (1, К) 
chord=xmax-xmin 
Gon? ЕСІ СЕ 
RT ЕКИ) спота 


Оооо ооо 


1=75 

j m 

k=km/2 

Vet ር. (cen Оаа ጋ 
gud kt gegen) 


ር: wIltei 7305) kk 25) 

do 20 k=1,km 
do. 10 IES 
i=ii(l) 
NS TOE 3070117 O 03 0 092 
нк, 

me 2211 =z (i, 7K) 

Ate (172300122122 (11) vel (11/velint,181.,5) 


ሮሆ እ ደደደ እደ ንን እእ ደቃን ንሙቷጥን እ ደን ሙንን አደ እቷልጴጊ፤።እ Ck Ck C CA CA Ck Ck Ok እሚ CK CK Ck ck Ck A Ck Ck ck ck Ck ok CA A Ck (Ck ck CA X 


ር. እእፈ ደፈ እ ፈጋ እጂ እእ እ እ እእ እእ እጂ እ እእ ፈጋ ጋፈ ቲን IK TI ቷን ን እእ እ AAA AAA AAA AAA 


КОО tormat(15,10f£10.:5) 


305 format ٦ И о 
610 formatí(' ***** warning *****',/, 
1 ANS RM 05. do mor 
mubelt./, 
2 Md kme eadi rom restart 112" 315) 
stop 
end 


Ó = es = مب‎ = = < <“ <“ =“ =“ - = =“ = - = = =“ = ) =“ = መ መ =“ ድሙ = መ መ =“ =“ =“ = መ መ =“ ሙመ =“ - - = =“ =“ =“ = < =< < < - ж- = نے‎ ዴዴ = - - = 


ӨӨ 


program shrink99 


This inputs data and multiplies it by 0.99 Тот а 7 
ЗИТ ЛК “factors 


> >x ው > 


real x(988),p(988) 


open (unit=2,file="'hotlv.dat'  statuss= 1 
open Junit#15,file=- "Iv.-tegın ር ር ንስ ከህ.) 
read(2,*) (x(1),1=1,988) 


do 20 i=1,988 
p(1)=x(1)*0.99 


- - = ወ a = - = ወ = <= - > = = - - = - — - = = - - - - - - ~ = = o u = - - - - = = = = = = = = > = = = = - - - - - - - 


20 continue 

А 

м -- print the data ont Ring) 

А 
Write рая | 
Write =m 
write (Шен ра 7067 
write boe 
Write ۴ )( 5 "9000ص‎ ٣ 
Write ls) 
write (15,2% ٦/٦ 
write» (15,95 ' 
write (15. 99 Ор лаје он о Боа 
write 115,%) 
write (1 መበ) (рај 71 
write (15) 
write (1534) (рауға БИ 
write (IS) 
write (ሞር ቁ (pi) 12512 564) 
write (L g" 
write (15, ۱۱-۹ 9ٍ 
write (157m 
Write (175, (pti) .1=618-67 0) 
Write (15,™) 
wrie (1579) (p (1); 1567215725) 
write (15 2777) 
write, (15%) о (а 1/2100) 
write (15 24 
Write (15у (Блу 1-2 Би 
wtite (157) 
write (15,*) ар) 
write (15,4) 
write (15,95. (PI) ICD 
write (15% 1) 
write 115,*) (p(1),1=936,988) 
end 
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